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Abstract 

We address two issues in the quantum electrodynamical description of photonic media with some disorder, neglecting material 
dispersion. When choosing a gauge in which the static potential vanishes, the normal modes of the medium with disorder satisfy a 
l" different transversality condition than the modes of the ideal medium. Our first result is an integral equation for optical modes such 
that all perturbation-theory solutions by construction satisfy the desired transversality condition. Secondly, when expanding the 
vector potential for the medium with disorder in terms of the normal modes of the ideal structure, we find the gauge transformation 
C~l that conveniently makes the static potential zero, thereby generalizing work by Glauber and Lewenstein [ Phys. Rev. A 43, 467 
(1991) ]. Our results are relevant for the quantum optics of disordered photonic crystals. 
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1. Introduction 

The quantum optics of random media is a young research 
field, studying the effects of randomness on quantum correla- 
tions and entanglement of quantum states of light in a multi- 
mode setting [ 1 , Traditionally, random media 
are studied with randomness against a homogeneous dielec- 
tric background. Recently researchers also realized that ev- 
ery real photonic structure, such as a photonic crystal, is in a 
sense a random medium, since there is inevitably some ran- 
domness on top of the ideal dielectric properties ||7l Ifa |9J] . The 
interplay between the randomness and an ordered inhomoge- 
neous dielectric background can sometimes be exploited. For 
example, in a photonic-crystal background slow light can pro- 
mote localization of light due to even minute random scatter- 
ing M m El Q . 

The typical starting point in the quantum optics of random 
media is the assumption of a multimode scattering matrix with 
elements subject to disorder HI 0, 51 . ^ ne l eve l °f mod- 
eling deeper is the quantum electrodynamics (QED) of these 
media, to derive the form of the scattering matrix and its de- 
pendence on the types of disorder in the medium. Here we 
aim to contribute to this QED description for spatially inho- 
mogeneous media with some additional disorder. For simplic- 
ity we assume that material dispersion can be neglected, as in 
Refs. Q d E EE HI, 

although more general quantized- 
field theories for dispersive and absorbing inho mog eneous di- 
electric media have also been developed lfl9l 20l 21, 22]. 

First we will derive a useful new integral equation for the nor- 
mal optical modes in a photonic medium with disorder. Related 
work on integral equations and Green-function methods can be 



found in Refs. EHIHI^I^IM^ and on dis- 

order in photonic media in Refs. jilliinisllMlillMIM^ 
4 lL 14211 . Instead of an integral equation for the modes involving 
a disorder potential and the usual Green tensor G of the unper- 
turbed medium, we introduce an alternative integral equation 
involving a kernel K 112711 that differs from G (details below). 
This K emerged naturally in a quantum optical description of 
light sources and scatterers in a photonic environment [27], and 
has since then been frequently employed in a quantum optics 
context, e.g. in Refs. II431 l44l l45l 14611 . Here instead we pro- 
pose a novel use to it, in an integral equation that we derive for 
optical modes of photonic media with disorder. We discuss its 
specific advantage that arbitrary-order perturbation-theory so- 
lutions automatically satisfy a desirable gauge condition. 

As our second topic we discuss an alternative to a normal- 
mode expansion, namely an expansion into modes that get cou- 
pled because of a perturbation. Disorder is such a perturbation. 
Several methods have been develo ped to describe disorder in 
photonic crystals l3lll3^lHl3lill4(Al4lll42ll . These methods are 
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useful both in the classical and in the quantum optics of disor- 
dered photonic crystals. Here our aim will be to generalize the 
coupled plane- wave description by Glauber & Lewenstein [15], 
to a coupled Bloch-mode description for example. Although 
there will be coupled modes in our theory, it is different from 
what is commonly known as a 'coupled-mode theory' 11471,1 1311 . 
which also finds application in quantum optics j48ll . 

The structure of this article is as follows: in Sec.|2]we briefly 
review the quantum electrodynamics of inhomogeneous dielec- 
tric media. In Sec. [3] a useful integral equation is derived for 
the independent optical modes, which is especially well suited 
for perturbation theory calculations of the effects of disorder on 
normal modes. Section |4] defines the problem to find a con- 
venient gauge transformation when starting from an expansion 
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into modes that are not the normal modes. This transforma- 
tion is constructed in Sec. [5] specified to the special case of a 
plane-wave expansion in Sec. [6] before we conclude in Sec. [7] 

2. Normal-mode expansion 

The quantum optical description of the electromagnetic field 
in a photonic medium with negligible material dispersion starts 
with the source-free Maxwell equations in matter V ■ D = 0, 
V-B = 0,D = VxH, and B = -V x E and the constitutive re- 
lations for a lossless nonmagnetic medium, D(r) = eoe(r)E(r) 
and B(r) = ^oH(r). We can introduce a vector potential A(r) 
and a scalar potential <5(r) such that 



the vector potential can be expressed as 



E = -VO - A, 
B = V x A, 



(1) 

(2) 



and where the dot denotes a time derivative. There is gauge 
freedom, i.e. one can make combined changes of the vector 
potential A — » A + and the static potential <D — > <t> - x that 
leave the electric and magnetic fields E and B unaltered. Here 
X(r, t) is an arbitrary scalar function of space and time. 

The usual steps from classical to quantum electrodynamics 
of photonic media are first to choose a convenient gauge, then 
to identify the canonical fields, next to express those fields into 
normal modes, and finally to associate non-commuting opera- 
tors with them 115, 3 18]. In the first step, the simplest gauge 
choice for lossless nondispersive photonic media is a general- 
ization of the free-space Coulomb gauge V ■ A = 0, namely the 
'generalized Coulomb gauge' 



V- [e(r)A] = 0. 



(3) 



We call such A 'generalized transverse' or 'e- trans verse', while 
Ref. J49ll uses 'quasi-transverse'. Advantage of this gauge 
choice is that the corresponding static potential O can be cho- 
sen identically zero, which leaves the vector potential as the 
only canonical variable in the Lagrangian density _£ that leads 
to Maxwell's equations, 



-C = ie £(r)[A(r, t)/c] 2 - [V x A(r, t)] 2 . 



(4) 



The field canonically conjugate to the vector potential equals 
eoe(r)A(r) = -D(r), and the Hamiltonian density becomes 



„ = M ' [VxA(r)f . 

2e e(r) 2/i 



(5) 



One can introduce commutation relations for A and -D directly, 
but a simpler equivalent procedure is to expand the fields into 
normal modes. The wave equation for the vector potential is 



_ _ . s(r)d 2 A n 
VxVxA+-V— = 0. 

c 1 ot 1 



(6) 



In terms of the complete set of normal modes f^(r) with mode 
index A that satisfy 




I [ a A (t)h(r) + ^(r)f» ] ■ (8) 



This is now a quantum mechanical operator where the creation 
and annihilation operators a [ and ax satisfy the usual harmonic- 
oscillator commutation relations. Apart from zero-point ener- 
gies that can be neglected here, the Hamiltonian is a sum of 
independent harmonic oscillators, 



H = 2 ftw 



(9) 



VxVxfir) + £ (r)^(r)=0, 

c L 



(7) 



This Hamiltonian defines the modes as 'normal modes' and 
leads to harmonic time dependence for the \ From the latter 
two equations and the relations (1112b . the mode expansions for 
the electric- and magnetic-field operators follow immediately. 



3. Gauge-respecting perturbation theory for normal modes 

In order to find the unknown normal modes of a medium II 
with dielectric function en(r), it is often useful to do this start- 
ing from the modes of another medium I, for which the nor- 
mal modes are either known or easier to compute or to inter- 
pret. Medium I is often an idealized structure with symmetries 
that make it easier to classify and find the normal modes, and 
medium II is its practical realization with some disorder. For 
example, all real photonic media have some disorder |Q,lil|9l]> 
unwanted or by design [50] or both ||l2ll . so that the realized 
dielectric function en(r) will be the sum of the dielectric func- 
tion of the ideal structure ei(r) plus an inevitable disorder term 
Ae(r) = fin(r) - £i(r). The normal modes {fn^} of the medium 
with disorder differ from the normal modes {fu} of the idealized 
structure, which for photonic crystals are Bloch modes. For the 
quantum electrodynamical description this does not pose any 
formal problems, for in principle one can follow the quantiza- 
tion procedure as discussed above, and write the vector poten- 
tial of the medium with disorder as 

A„(r, t) = £ [ Mfl«Wr) + a L( f )fiu(r) ] , (10) 

The modes {fn^} and the vector potential An satisfy the gauge 
condition V • [ £n(r)Arj(r) ] = 0, and this is the gauge in which 
the static potential for medium II vanishes. 

In practice, it can be quite a challenge to find the normal 
modes of the vector potential for a complex photonic medium 
with disorder. Here we develop a calculational tool that can 
make it simpler. We consider the problem how to find a normal- 
mode expansion for the vector potential of medium II, given the 
normal-mode expansion of Ai of medium I in terms of mode 
functions fu with V ■ [ei(r)fu(r)] = 0. For medium II we also 
wish to work in the analogous convenient gauge in which V • 
[en(r)An(r)] = so that the static potential vanishes. Our goal 
is now to find the normal modes fm. We show that the standard 
perturbation expansion does not have the desired transversality 
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property, and propose a new type of perturbation expansion for 
the normal modes of the vector potential that resolves this issue. 

First we derive in a few lines the standard integral equation 
for the modes of medium II. There is a trivial way to rewrite 
their defining equation (0 as 



- V x V x %(r) + ei(r)-f f u „(r) = V(r, a, J ■ f^(r), (11) 

with the idealized dielectric function ei(r) on the left-hand side, 
and on the right the term with the perturbation potential 

V(r, to) = -[en(r) - £i(r)](w/c) 2 l = -Ae(r)(co/cf\, (12) 

where I is the unit tensor. Next we introduce the Green ten- 
sor Gi, defined as the solution of the wave equation with delta- 
function source term, 



-VxVxG 1 (r,r»+ £l (r)-iGi(r,r» = S(r-r% (13) 

c z 

There is an analogous equation for Gn. By combining the last 
four equations, we can find a mode fm starting with an unper- 
turbed mode fu, using the exact Lippmann-Schwinger integral 
equation 

fiu(r) = fu(r) + j dr' G^r, r', co A ) • V(r', co A ) • f m (r').(14) 

The fiju have the same label A, because the two are related by 
this integral equation. 

We know that V ■ [£i(r)f u (r)] = and from Eq. CCD that 

V ■ [eii(r)fm(r)] should vanish for the normal modes of medium 
II, i.e. for the solution fyy of the integral equation ( fT4t . 
However, in zero-order perturbation theory one has the solu- 
tion f^j = fu that clearly breaks the gauge condition of £rr 
transversality. Likewise, in first-order perturbation theory or 
Born approximation, where one replaces fru by fy within the 
integral in Eq. ( fT4t . one finds an improved approximation 
that nevertheless breaks the gauge condition for fru- The con- 
dition is satisfied by the exact (infinite-order) solution only. 
Now doing perturbation theory up to some finite order of course 
means that some controlled error is introduced. However, an al- 
ternative integral equation that generates solutions that satisfy 

V • [en(r)ffu (r)] = for all orders of approximation n is clearly 
to be preferred. Below we present just that. 

The Green tensor Gi of Eq. ( fT3l can be written as the sum of a 
generalized transverse part G^ and a longitudinal {i.e. curl-free) 
part G\, which can be expanded in a complete set of generalized 
transverse modes of medium I as 112711 



GT(r,r» = c 2 J] 
G^(r,r',w) 



f u (r)f^(r') 



(15) 



A (M> + i^) 2 - w 2 ' 

Si(r)(co/c) 2 (cj/c) 2 

From this it is easy to see that the total Green function can al- 
ternatively be expressed as 112711 



where we introduced the tensor Ki that has the mode expansion 
_\2 f u (r) f u (r') 

\ ti) I -4- ir 



is , ' \ 2^ tu(T 



{to + irf) 2 - toj 



(18) 



An important difference between K and G is that V r • 
[si(r)Ki(r, r', to)] = 0, as follows from the mode expansion 
Eq. (flST l. whereas V r • [£i(r)Gi(r, r', to)] ± 0. This same Green 
tensor K emerged naturally in the multiple-scattering formalism 
of light interacting with atoms in photonic media in Ref. [27], 
where there was no way around it, so to say. Here instead we 
just choose to rewrite the integral equation (TBI in terms of the 
new Green tensor, because it produces a more convenient in- 
tegral equation. After using the identity Eq. dT7T > to make the 
replacement of G by K in the integral equation (TT4l > and evalu- 
ating the integral over the delta-function term of Eq. ( fTTl i with 
the help of the expression (fT2l for V, we obtain after rearrang- 
ing the exact integral equation 

enWfmW = fii(r)f u (r) (19) 

+ Sl (r) J dr' Kj(r, r', cod ■ V(r\ to*) ■ f IL1 (r'). 

Both terms on the right-hand side are divergence-free, the first 
one by assumption, and for the second term it follows from the 
ei-transversality of Ki as discussed above. It follows that the 
left-hand side is also divergence-free, so that fru( r ) on the left 
is indeed a generalized transverse normal mode for medium II. 
Even better, all finite-order perturbation-theory solutions based 
on this integral equation ( fT9l l have this same property. This is 
obvious for the zero-order solution fr??( r ) = £i(r)fu( r )/ £ ii( r X 
found by putting V to zero in Eq. (fT9l . And to first order in the 
perturbation potential V, i.e. in Born approximation, we have 



£n(r)0r) 



£i(r)f M (r) (20) 
£l (r) J dr' K^r'.wJ-Vfr'.^'Or'). 



Gi(r,r',w) = Ki(r,r',a») + 



6(r - r')l 

e!(r)(w/c) 2 ' 



(17) 



The left-hand side of this equation is divergence-free, for the 
same reasons as given above for the left-hand side of Eq. dT9b 
for the exact (but implicit) solution. 

The new integral equation ( fT9l is advantageous as compared 
to Eq. (TBI for another reason. Often one would like to write the 
solution as a linear combination of solutions of the unperturbed 
system. In Eq. dT4l > one would be tempted to write the as lin- 
ear combinations of the fi v , but such an expansion would be in- 
complete because of the different transversality relations of the 
two types of modes. In Eq. ( fT9l one can define bqju = eijifyu, 
which are divergence-free (just like the transverse plane waves 
in which they could be expanded). The hiu can therefore be 
completely expanded in terms of the hy, or hiu = Zv Ci v hi v , 
and hence the fiu in terms of the fu, ei(r), and en(r). 

To give a simple example of the advantage of our expan- 
sion, assume that a smooth background ID dielectric function 
£b(z) is perturbed and in a finite z-interval is replaced by the 
smooth dielectric function sa(z), with discontinuities on the in- 
terfaces. At an interface, it is well known that the tangential 
components of the electric field and the normal components of 
the displacement field are conserved. However, already if we 



use the zero-order solution of Eq. (fT4l . the tangential compo- 
nents of the electric field are correctly continuous, but the nor- 
mal components of the displacement field are not. Conversely, 
the zero-order solution Eq. d20b incorrectly gives discontinuous 
parallel components of the E-field but correctly gives continu- 
ity of the normal component of the D field. So neither zero- 
order solution is the exact solution of course, but the advantage 
of the zero-order solution Eq. ( l20t is that it correctly satisfies 
the gauge conditions V • [eb(z) A(r)] = in the B-region and 
V • [sa(z)A(t)] — in the A-region. In this first example we 
elaborated on the lowest-order approximation. However, the 
advantage in quantum electrodynamics of our gauge-respecting 
series ( fT9l ) exists for every order of perturbation theory. 

It is also instructive to do an exact calculation with the 
new integral equation ( fT9l ), of local-field effects for example. 
The electric field in an infinitely small spherical empty cavity 
around a point-like emitter at r = embedded in a homoge- 
neous medium with dielectric function e > 1 is enhanced by a 
so-called local-field factor L, see 115 ll kill and the derivation 
in IU5I1 based on Eq. ( fT4b . so that the spontaneous-emission 
rate becomes L 2 \/e times the free-space value. Assumed is 
that the field inside the sphere is constant [f(r) = f(0)] and 
important to know is that only the delta-function contribution 
6(r - r')l/(3e( oj/c) 2 ) of G contributes to the volume integral in 
the small sphere, see 15211 . Eq. (fT4l then immediately gives L = 
1 + (e— l)L/(3e), which has as solution the well-known empty- 
cavity local-field factor L = 3s /{2s + 1). Likewise one can 
derive a local-field factor from the new equation ( TT9l >, knowing 
from Eq. (fTTT i that the delta-function contribution of K(r, r', u>) 
for the homogeneous medium is -25(r-r')l/(3e(w/c) 2 ): hence 
the integral equation ( fT9T > gives L-s - s(s- l)2L/(3e). Again 
we find the L = 3s/ (2s + 1 ). This illustrates that the two pertur- 
bation series are equivalent, in the sense of producing identical 
exact (infinite-order) solutions. Nevertheless, finite-order ap- 
proximations of the equivalent series will generally differ. 

4. Gauge problem for non-normal mode expansions 

For the QED of photonic media with some disorder, either 
one chooses to expand the field operators in terms of the normal 
modes {fn^,}, as above, or one characterizes the disorder by the 
way the idealized modes {f^} become coupled. 

We will now turn to the latter approach and expand the vector 
potential An of medium II in terms of the normal modes {fu} of 
medium I. This is certainly possible, but then An satisfies the 
gauge condition V-[ £i(r)An(r) ] = 0. With that gauge condition 
for an en medium, a nonzero static potential is required to make 
the displacement field divergence-free. Since our goal is to set 
up a canonical theory in terms of only the vector potentials, a 
gauge transformation is required. Glauber & Lewenstein lfl5ll 
pointed out this complication of non-normal mode expansions 
in QED when considering the special case ei(r) = 1 . 

We therefore ask whether a gauge transformation exist 
that transforms vector and static potentials that satisfy V ■ 
[ £i(r)A(r) ] = and O + into vector and static potentials 
satisfying V • [ en(r)A(r) ] = and <I> = 0, respectively. Glauber 
& Lewenstein used properties of the free-space transverse delta 



function to find the required gauge transformation for the spe- 
cial case £i(r) = 1 [15]. With applications to disordered pho- 
tonic crystals in mind, we will now construct a more general 
gauge transformation for arbitrary smooth positive dielectric 
functions £i,n(r). We thereby generalize Ref. IU5I1 on this point, 
and our method will be different. 

As an Ansatz, we consider gauge transformations of the form 



(21) 



X(r, t) = V Jt— [flu(fki(r) + h.c. 



We wish to show that a particular ^(r, f) exists such that the 
vector potential of the en-medium can be expressed as 

A u (r,f) = J] ^^^{ «u(0[fu(r)+ Y^O-)] +h.c. }, (22) 

satisfying the generalized Coulomb gauge condition V ■ 
[fin(r)An(r, f)] = 0. This gives for every separate mode A 



V ■ {enW [fu(r) + V^(r)] }= 0. 



(23) 



This gauge condition defines the 'gauge problem' for non- 
normal modes, and we solve it in the next section. 



5. Gauge transformation for non-normal modes 

In Sec. [3] we expressed the modes of the realized structure 
II in terms of the modes of the ideal structure I. For our pur- 
pose of finding the gauge transformation x, we will consider 
the less commonly used perturbation expansion in the other di- 
rection, expressing the modes of the ideal structure in terms of 
the modes of the realized structure. The wave equations for the 
ideal modes fi can be written as 

- V x V x f u (r) + e n (r)^f u (r) = -V(r, oj a ) • f u (r), (24) 

and notice the minus sign on the right-hand side, as compared 
to Eq. (fTTT i. Again an implicit solution in terms of an integral 
equation can be given, 

f u (r) = f M (r) - J dr' G u (r, r', co A ) ■ V(r', co A ) ■ f u (r').(25) 

This exact relation is interesting for our purposes, since we are 
looking for a gauge transformation of the left-hand side of this 
equation that makes it en-transverse, and the first term on the 
right-hand side of this equation already has this property. The 
last term does not, which is consistent with the left-hand side 
being ei-transverse. 

At this point we use a property of the Green tensor Gn that we 
call its generalized Helmholtz decomposition: it has a unique 
decomposition into an en -transverse part Gj and a longitudinal 
part G n , or G u = G£ + GJ}, with 

V • [e u (r)G„(r, r')] = 0, and V x G^(r, r') = 0. (26) 



The proof of this property of Green tensors is given in Ref. [27], 
and is based on a similar unique decomposition of vector 
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fields IU7L [18Q - We can make direct use of this Green-tensor 
decomposition, simply by adding a term J G}j ■ V ■ on both 
sides of Eq. (|25T >. In the same short-hand notation, we obtain 



(27) 



term d28l l and the expansion (fT6b of the longitudinal Green ten- 
sor Gn into the normal modes fnko-, we find 



fu + f(3k ■ V ■ i u = f lu - jGl ■ V ■ f u , 



It is now obvious that the right-hand side is en-transverse, so 
we conclude the same for the left-hand side. Moreover, the sec- 
ond term on the left is longitudinal. Thus we find what we set 
out to prove, namely that for every mode function fu, there ex- 
ists a gauge transformational that satisfies the gauge condition 
Eq. d23l . But we know more than its mere existence, for we 
find that the gauge transformation has the particular form 



V^(r) 



(r,r',^)-V(r', W/1 )-f u (r'), 



(28) 



where G n can be expressed in terms of the modes fru, analo- 
gous to Eq. dl6l . 

It is useful to summarize what we have achieved here. Start- 
ing with a fi-mode expansion of the vector potential (which 
thereby was ei-transverse), we have found the gauge transfor- 
mation (l28T i that makes it en-transverse. This is the gauge for 
which the static potential of the en-medium can be chosen iden- 
tically zero. Since after the gauge transformation the vector po- 
tential is the only canonical field, it is a simple matter to express 
the electric and magnetic fields, using Eqs. ( 1 1 12b : 



E u (r,r) = "Z J^(Wf)[fu(r) + V^(r)]+h.^ 

Bu(r '° = 2 ^2^y auit)[V ' xf u( r )] +h - c -}> 00) 

where the gauge term in the magnetic field vanished, being the 
curl of a gradient. The fact that in this gauge the magnetic 
field of medium II is expanded in the normal mode functions 
of medium I can be numerically advantageous. The electric 
field has an additional non-vanishing gauge term, taking over 
the role of - and mathematically identical to - the static poten- 
tial in the initial gauge. The displacement field Dn is given by 
eoen(r)En(r), and hence is divergence-free as it should. 

6. Special case: plane waves as non-normal modes 

Here we study how our results of the previous section sim- 
plify in the special case that ei(r) = 1, in other words if for 
the ideal disorder-free medium I we choose free space, with its 
transverse plane-wave modes exp(ik • r)e kcr , where <x = 1,2 
labels the two orthogonal unit vectors e kcr perpendicular to the 
wave vector k. In medium II with dielectric function en(r), 
these transverse plane waves are not the normal modes of the 
vector potential An- But we can just expand An in transverse 
plane waves, and use the gauge transformation^ that we found 
in the previous section to end up in the gauge in which the 
static potential vanishes. Using the explicit form of the gauge 



Y*w(r) = 



1 - ea(r) 



£n(r) 



e lkr e kff 



(31) 



+ J] fm«r(r) f dr' f IIk<r (r') • [1 - e u (r')] e lk r 'e k(T . 
ko- J 

Therefore, the spatial dependence [fu(r) + Vy,i(r)] of the vector 
potential An in Eq. ( 122b in our special case becomes 

— ^ + 7 !Wr) f dr'Wr>[l-e„(rO]e ikr 'e kcr .(32) 

It follows that the vector potential indeed satisfies V • 
[en(r)An(r)] = 0, that the electric field En = -An satisfies the 
same condition, and that Dn = eoen(r)En is divergence-free. 

For the magnetic field Bn it follows from Eq. (130b that it can 
be fully expanded in the free-space transverse plane waves. The 
displacement field has the same spatial dependence as An in 
Eq. (l32l . but multiplied by en(r). The result is a transverse 
plane wave, plus the divergence-free gauge term that is not a 
plane wave. By Fourier-expanding the gauge term into trans- 
verse plane waves and regrouping, one could alternatively ex- 
pand the displacement field into transverse plane-wave modes, 
which is the approach of Glauber & Lewenstein [15]. 



7. Conclusions 

After a short review of the quantization of the electromag- 
netic field in inhomogeneous dispersionless dielectrics, we ar- 
gued that a convenient choice of gauge is the one for which the 
static potential vanishes, leaving the vector potential as the only 
dynamical field. In such a gauge, a mode function of the vector 
potential for an idealized structure has different transversality 
properties than for a realistic structure, due to disorder. 

Perturbative solutions based on a standard integral equation 
for modes in disordered photonic media do not give the desired 
transversality property. Here we proposed an improved integral 
equation, in terms of the Green tensor K [27] rather than G, that 
does give such solutions to any order of perturbation. This can 
become a useful numerical tool to find exact or approximate 
normal modes, for example in photonic crystals with disorder. 

Besides normal-mode expansions, one can alternatively ex- 
press the field operators in terms of modes that are not indepen- 
dent, with couplings amongst them due to the disorder. Start- 
ing with an expansion of the vector potential of a disordered 
medium in terms of modes of an idealized structure, we con- 
structed the gauge transformation that makes the static potential 
identically zero, and obtained expansions for the field operators 
into modes of the idealized structure plus gauge terms. 

We focused on issues related to the choice of a gauge in the 
quantum electrodynamics of photonic media. Our results will 
be useful for developing a full QED theory of modes in complex 
photonic systems that become coupled due to disorder. 
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